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Abstract 

Recently, we have introduced a unified theory to deal with perturbative and 
non-perturbative QCD contributions to hadronic structure functions in deep 
inelastic scattering. This formulation is realized by combining the coordi- 
nate space approach based on light-front current algebra techniques and the 
momentum space approach based on Fock space expansion methods in the 
Hamiltonian formalism of light-front field theory. In this work we show how 
a perturbative analysis in the light-front Hamiltonian formalism leads to the 
factorization scheme we have proposed recently. The analysis also shows that 
the scaling violations due to perturbative QCD corrections can be rather easily 
addressed in this framework by simply replacing the hadron target by dressed 
parton target and then carrying out a systematic expansion in the coupling 
constant Og based on the perturbative QCD expansion of the dressed par- 
ton target. The tools employed for this calculation are those available from 
light-front old-fashioned perturbation theory. We present the complete set of 
calculations of unpolarized and polarized deep inelastic structure functions 
to order a^. We extract the relevant splitting functions in all the cases. We 
explicitly verify all the sum rules to order Og- We demonstrate the validity of 
approximations made in the derivation of the new factorization scheme. This 
is achieved with the help of detailed calculations of the evolution of structure 
function of a composite system carried out using multi-parton wavef unctions. 



PACS numbers: 12.38.-t, 13.85.Hd, 13.88.+e, ll.30.Rd 
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Typeset using REVT^ 



I. INTRODUCTION 



A general cross section in hadron physics contains both short distance and long distance 
behavior and hence is not accessible to perturbative QCD. Factorization theorems allow 
one to separate the two behaviors in a systematic fashion. Usually the short distance (per- 
turbative) properties are calculated with Feynman diagrams where the most popular choice 
of regulator is dimensional regularization. So far there is no non-perturbative implementa- 
tion of dimensional regularization in field theory. The long distance (non-perturbative) part 
is given in terms of a set of operator matrix elements which are left for computation, for 
example, in lattice gauge theory. Since factorization scale is an artifact reflecting our present 
inability to do computations in QCD, the two sectors (perturbative and non-perturbative) 
should merge smoothly. Since currently available formalisms employed to tackle the two 
sectors use different regulators, different degrees of freedom, etc., this goal is difficult to ac- 
complish in practice. It is desirable to have a method of calculation where same formalism 
is used to deal with both perturbative and non-perturbative regions of QCD. 

Recently we have proposed a new method of calculation of deep inelastic structure 
functions that combine current algebra techniques and Fock space expansion methods in 
light-front field theory in a Hamiltonian framework. We have arrived at expressions for 
various structure functions as the Fourier transform of hadron matrix elements of different 
components of bilocal vector and axial vector currents on the light-front. By expanding 
the state of the hadron in terms of multi-parton wave functions, non-perturbative QCD 
dynamics underlying the structure functions can be explored. We have also presented a 
novel analysis of power corrections based on light-front power counting. 

In this work we show that a perturbative analysis in the light-front Hamiltonian formal- 
ism leads to the factorization scheme proposed in Ref. 0. The analysis also shows that 
the scaling violations due to perturbative QCD corrections can be rather easily addressed 
in this framework by simply replacing the hadron target by dressed parton target and then 
carrying out a systematic expansion in the coupling constant as based on the perturbative 
QCD expansion of the dressed parton target. The tools used are those belonging to light- 
front old fashioned perturbation theory which employs transition matrix elements and 
energy denominators instead of Feynman propagators. 

The advantage of using the light-front Hamiltonian formulation is that the matrix ele- 
ments can be naturally defined in the light-front gauge (^4+ = in light-front coordinates). 
With this gauge choice, there is no need for the path ordered exponential between fermion 
field operators in the bilocal current which is mandatory in other popular gauge choices. 
Since we do not have to deal with four dimensional integrals involving Feynman propaga- 
tors we do not encounter some of the problems associated with the use of the non-covariant 
gauge condition = in usual covariant perturbative QCD calculations. 

Meanwhile, the evaluation of the matrix elements in our approach (see sections IV and 
V) is straightforward, which also greatly clarify the physical picture of DIS. For example, 
matrix elements of the transverse component of the bilocal vector and axial vector currents 
are easy to analyze in the present method. In sharp contrast, there are well-known problems 
IP associated with 75 in dimensional regularization. Also presence of quark masses poses 
no problem in our calculations which has been bothersome for a long time in the standard 
OPE or Feynman diagram approach of DIS. 
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Since our approach deals with probabihty amphtudes, real and virtual processes are 
calculated to the same order without any difficulties, in contrast with the Altarelli-Parisi 
method which deals with probability densities. We ffist present result for the unpolarized 
structure functions F2{x) for dressed quark (extracted from the plus component of the 
vector bilocal) and gluon target to order as- We extract the relevant splitting functions. In 
these cases we also explicitly verify the longitudinal momentum sum rule. We also present 
result for the polarized structure function gi for a dressed quark. Furthermore, interference 
effects are straightforward to handle. As a result, we explicitly show the invalidity of the 
popular twist classification [Q. We demonstrate this by analyzing the matrix element of the 
transverse component of the bilocal vector and axial vector currents (see Sees. IV and V 
for details). In the conventional OPE analysis P], it is customary to ignore the non-trivial 
structure of the state, and consider the operator structure alone to draw conclusions. A 
case in point is the transverse component of the bilocal vector current. This operator is 
not diagonal and in the twist analysis of Ref. [^] will appear to have twist three. Hence it 
appears, as if this operator matrix element cannot have parton interpretation. However our 
explicit calculations show that this operator matrix matrix element is indeed twist two and 
has the familiar parton interpretation. This becomes clear only after the evaluation of the 
matrix elements which includes off-diagonal ones [ITU . 



For a second example, in the case of the transverse polarized structure function, it is 
popular to ignore quark mass and stress quark-gluon correlations. We show by explicit 
calculations that the operators that involve 75 which do not have explicit dependence on 
the quark mass ruq, come out proportional to rriq when matrix element is taken between 
dressed quark states. In particular all the operators contribute at the same level to the 



evolution of gx- This is shown to be essential |TT[ for the g2 structure function to obey 
Burkhardt-Cottingham sum rule []T2| in perturbative QCD. 

In this work we also provide a detailed calculation of the evolution of F2 structure function 
for a hadronic bound state. Entire analysis is carried out using multi-parton wave functions 
in momentum space. Both real and virtual processes are accounted for and in the lowest 
order analysis one begins to see the emergence of the standard evolution equation. The 
detailed analysis justifies the approximations made in Sec. II which lead to factorization 
to all orders in perturbation theory. Then what remains unsolved is the nonperturbative 
contribution that is defined in the same framework. At present, the nonperturbative QCD 
dynamics and hadronic bound states on the light-front are indeed explored in this same 
framework. Therefore, we hope that the present work (ref. |0] and the present paper) can 
really provide a natural connection of the fully theoretical understanding to the experimental 
phenomena of DIS. Note that for exclusive processes, factorization has been proven using 
light-front formalism by Brodsky and Lepage ||^. 

The plan of this paper is as follows. In section II we present a perturbative analysis in 
the light-front Hamiltonian formalism which leads to the factorization and to the concept 
of the structure function of a dressed parton in DIS. The tools necessary for the evalua- 
tion of these functions are discussed in Sec. HI. Unpolarized and polarized dressed parton 
distribution functions are discussed in Sees. IV and V respectively. We also present the 
explicit verification of the appropriate sum rules in these sections. A detailed analysis of 
the structure function of bound states is carried out in Sec. VI which justifies a posteriori 
the approximations made in the study in Sec. II. Finally Sec. VII contains discussion and 
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conclusions. 



II. FACTORIZATION: A PERTURBATIVE ANALYSIS 

In this section we show in detail, how the factorization picture discussed in Ref. 
emerges in a perturbative analysis carried over to all orders in the case where the bilocal 
operator involved does not change the particle number. The analysis leads to the concept 
of the structure function of a dressed parton in DIS. 

To explicitly demonstrate the factorization picture on the light-front, we consider the F2 
structure function as a specific example in this section. For simplicity we drop reference to 
the flavor, then 

^^^^^= ^ / di-e-'^^^^-^{PS I [(^+)t(r)^+(0) - (^+)t(0)^+(r)] I PS). 

(2.1) 

From the discussion in Sec. V.D of 0, we have 



X 



q{x,Q') = —J dre-^^^«"^ E (^V I n,){n2 \ PSfi') 



{m I f/,-^[(^+)t(r)^+(0) - (^+)t(0)^+(r)Jf/h I n^) , (2.2) 

where Uh = T+exp | — | J-00 '^^^-Pint(^^)}) ^'^d Pj^^ = P~^^ + P~^f^ is denoted as the hard 
and mixed light-front interaction Hamiltonian. 

In the following, we shall not explicitly evaluate contributions from intermediate states 
which involve vanishing energy denominators. These contributions are most conveniently 
included by introducing the wave function renormalization constant associated with the 
parton active in the high energy process. 

Let us consider first few terms in the expression given in Eq. ( p.2|) . The lowest order 
term yields the function 

q^'\x,Q')=q{x,^,')=Y.\{PS,^^Un)? 

n 

= Y: rd'k^iPSfi' I bl{yP^,k^MyP^,k^) I PV) • (2.3) 

s •' 

Terms linear in the interaction Hamiltonian vanishes since the plus component of the bilocal 
operator conserves particle number on the light-front. Consider the second order contribu- 
tion. 



q^^\x,Q')= ^ (^V I n) f dxt{n \ p-,(x+) | m)(m | O \ p) 

(p| f dxtPUxt)\k){k\PSfi'). (2.4) 



Here we have denoted 
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(^ = ^1 cire-^^^^"^[(^+)^(r)^+(o) - (^+)^(o)^+(r)] (2.5) 



Using, 



we have, 



P^nt{^^) = ei''Fr.M)p-^(o)e-i^^.M)^ (2.6) 



{m\0\p){p\ P,r,{0) I A;)(PV I I PV)- (2-7) 

The states | n), and | /c) are forced to be low energy states with (/c"*")^ < /x^. We can 
restrict the states | m), | p) to be high energy states with (/c"*")^ > jj,"^. The bilocal operator 
O picks an active parton in a high energy state whose longitudinal momentum is forced to 
be xP^. Further we need to keep only terms in Pj7it which cause transitions involving the 
active parton. (Transitions involving spectators lead to wave function renormalization of 
spectator states which are canceled by the renormalization process as shown explicitly in 
Sec. VIB.) 

To order ctg, a straightforward evaluation (see later) leads to 

q{x, Q') = M {q{x, I,') + ^Cfln% f ^P(x/y)q{y, fx')} (2.8) 

where Af is the wave function renormalization constant of the active parton and P is the 
splitting function. Including the contribution from the wave function renormalization con- 
stant to the same order (ccs), we get, 

q{x, Q2) = jdy V{x, Q2; y, i/) q{y, /x^), (2.9) 

where the hard scattering coefficient 

V{x, = 5{x - y) + ^Cfln^ / dz6{zy - x)P{z) (2.10) 

ZTT IX'^ Jo 



with P{x) = P{x) - 5{l - x) /J dyP{y). 

We note that the above analysis can be carried over to all orders in perturbation theory 
with the result 

V{x, Q^- 2/, 11")= {yP+, k^, s I Ui;'OU, \ yP+, k^, s) 

= {yP^, k^, s; (dressed) \ O \ yP+, k^, s; (dressed)), (2.11) 

In evaluating the above expression, only in the interaction Hamiltonians in the extreme 
left and extreme right of the time ordered product we need to keep mixture of soft and hard 
partons. This is governed by P^^^ ■ They are needed to cause of transition of a soft parton 
to a hard parton. In the rest of the interaction Hamiltonians occurring in the chain, the 
partons are restricted to be hard, i.e., they are determined by Pj^^ only. For the leading 
logarithmic evolution we are discussing here, they appear ordered in transverse momentum. 
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III. TOOLS OF CALCULATION 



In this section, we outline the basic tools for calculating the perturbative contribution 
to the structure functions, namely the hard scattering coefficients, V{x,Q^]y,fi'^), given by 
Eq. ( |2.11|) . If we set k = P, then y = 1 and the hard scattering coefficients just become the 
structure functions of dressed quark and gluon targets in DIS, 

/f (x, Q') = ^ / dve-'"^' p{ks\ [V^(r)r.^(0) T h.c] \ks), . (3.1) 

As a matter of fact, we can compute the perturbative QCD correction to the hadronic 
structure functions by calculating the structure functions of the dressed quarks and gluons. 
In old-fashioned light-front perturbation theory, the dressed quark or gluon states can be 
expanded as follows @]: 

I Ps),= U,\Ps) = VJ7{ I Ps) + ± I 

+ y\Tn)- '-^^ '-^ vn^_l_ — ^ + ..1 (3.2) 

' P- - P~) P- - P-) J ^ ^ 

mn V-' ^ mJ\ n J 

where |Ps), the bare single particle state, and | n), the two-particle state, | m), the three- 
particle state, etc., are eigenstates of the free Hamiltonian. Single particle state \Ps) is 
omitted in the sum over the states in the above expansion. Introducing the multi-parton 
amplitudes (wave functions), 

(n I Hr^ I Ps) 
" (P-_P-) ' 

^ _^ {m\p-^\n){n\Pr^\Ps) 

the expansion in Eq. (|3.2| ) takes the form 

I Ps)^ = VU[ I Ps) + ^ $ J n) + ^ I m) + ...} . (3.4) 

n m 

In the above expressions, Pj^f^ and P~^^ are the interaction parts of the canonical light- 
front QCD Hamiltonian, but the former contains the mixed soft and hard partons and the 
latter only has hard partons. The canonical light-front QCD Hamiltonian is given by the 
following form in our two-component formalism [Q: 

^int = J dx d X ^T^qqg + 'Hggg + 'Hqqgg + T^qqqq + ^^gggg^ ) (3-5) 



and 



[a ■ d'^ + m) 



+ [-^){a-d^-m)a-A^y, (3.6) 
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999 



u 



d+ 



9^ 



9^ 



7^, 



2^? 



9+ 



d+ 



H 



9999 



^^9999^ "I" "^gggg"^ ' 
where the dynamic fermion field operator 



with 



X 2(27r)VA;+ 
and the transverse gluon field operator 

dk+cPk 

X 

with 







bxik)e-''''' + dl^{k)e'''^ 



{6a(A;),6UA;')} = {dx{k),d{,k'] = 2{2T^fk+5{k+ - k'^)5''{k^ - k' 



ax{k), a{,{k') = 2{2nf k+6{k+ - k'^)5\k^ - k^') 



(3.7) 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



(3.12) 



(3.13) 



(3.14) 
(3.15) 



and xa is the eigenstate of in the two-component spinor of ipj^ by the use of the following 
light-front 7 matrix representation [|TB| , 
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"0 -i" 




"0 i' 




' -ia^ 


i 


, 7^ = 


i 


za' 



(3.16) 



with = a^, (T^ = — cr^) and £*(A) the polarization vector of transverse gauge field. 



IV. UNPOLARIZED DRESSED PARTON STRUCTURE FUNCTIONS 

In this section we present the calculations of the F2 structure function for dressed quark 
and gluon targets. 
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A. Dressed quark structure function from the plus component 



The F2 structure function of a dressed quark is given by 



X 



1 



4nP+ 

The gluon structure function fl]] is defined by 

1 



47rP^ 



dr]e-'''\{k\ I {-)F+'"'ir)F^t{0) I kX)p. 



(4.1) 



(4.2) 



The dressed quark or gluon state can be obtained by the perturbative expansion in the 
old fashioned time-ordered Hamiltonian formulation, as given by Eq. ( |3.2| ). But we can also 
find such states by solving the light-front bound state equation. Let us take the state | P) 
to be a dressed quark which obeys the eigen value equation 



i=l 







^gg 





{q\H^nM) {(l\Hgt\q9) 



^g 



gg 



(4.3) 



Explicitly, expanding the state in terms of bare states of quark, quark plus gluon, quark 
plus two gluons, etc, we have. 



Pa), 



fM,\h\P,a)\Q) 



+ E 



dk2 d ko I / T-> I J 7 \ \ 
--i'2{P,o- I fci,Cri;fc2,A2j 



X J2{27,y^P+6'{P -ki- k2)¥{ki, ai)a\k2, X2) \ 0) 



(4.4) 



The factor J\fq is the wave function renormalization constant for the quark and the function 
ip2iP,o' I kiai, k2X2) is the probability amplitude to find a bare quark with momentum ki 
and helicity ui and a bare gluon with momentum k2 and helicity A2 in the dressed quark. 



Introduce the Jacobi momenta (xj, nf- 



so that 



(4.5) 



(4.6) 



The amplitude t/'2 is related to the amplitude $2 in Eq. 

y^i^2{kt,ki)=^2{x,,^^i 



by , 



(4.7) 
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The two particle amplitude V'2 is given by 



i{j2{P,(^ I Pu si;P2,p2) = 
9 1 



P+ vf 



X 



2(27r) 



rpa 



+ + 

P2 Pi 



im 



a^.P-^ + im 



-a 



a 



P2' 



(4.8) 



We rewrite the above equation in terms of Jacobi momenta {pf = XiP^, i^i = pi + XiP-^) 
and the wave functions $j which are functions of Jacobi momenta. Using the notation 
X = Xi, Ki = K, and using the facts Xi + X2 = 1, ki + k,2 = 0, we have 



n,P2 



(^X ^ /tj j X X^ Aij ) 



1 



^2 _ mf+Olf)^ 



X 



2(27r) 



- 2- 



K 



1-x 

a-^.K,-^ — im 



-a 



a im 



X 



X 



P2' 



(4.9) 



Evaluating the expression in Eq. ( |4.1| ) explicitly, noting that in the present case the 
contribution from the second term in this expression is zero, we get the quark structure 
function of the dressed quark 



X 



— X 



+ ^ y dx2 j dPl^i j (i^K^(5(l — X — X2) 



(4.10) 



This equation makes manifest the parton interpretation of the quark distribution function, 
namely, the quark distribution function of a dressed quark is the incoherent sum of prob- 
abilities to find a bare parton (quark) with longitudinal momentum fraction x in various 
multi-particle Fock states of the dressed quark. Since we have computed the distribution 
function in field theory, there are also significant differences from the traditional parton 



model |T^. Most important difference is the fact that the partons in field theory have trans- 
verse momenta ranging from zero to infinity. Whether the structure function scales or not 
now depends on the ultraviolet behavior of the multi-parton wave functions. By analyzing 
various interactions, one easily finds that in super renormalizable interactions, the transverse 
momentum integrals converge in the ultraviolet and the structure function scales, whereas 
in renormalizable interactions, the transverse momentum integrals diverge in the ultraviolet 
which in turn leads to scaling violations in the structure function. 
Taking the bare and dressed quarks to be massless, we arrive at 



L\ |2 



{2tt 



1 — X 



(4.11) 
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where Cf = ^2n^ ■ Recalling that | $2(x, k^) p is the probability density to find a quark 
with momentum fraction x and relative transverse momentum k-^ in a parent quark, we 
define the probability density to find a quark with momentum fraction x inside a parent 
quark as the splitting function 



c 



l + x'^ 



f- 



X 



(4.12) 



Clearly, the probability density to find a gluon with momentum fraction x inside a parent 
quark is defined as the splitting function 



Gq 



x] 



C-. 



l + (l-x) 



X 



(4.13) 



The transverse momentum integral in Eq. ( [4.11[ ) is divergent at both limits of integration. 
We regulate the lower limit by /i and the upper limit by Q. Thus we have 



X 



In —r 



x) + —Cf 

' 27: ^ 1-x 



1^- 



The normalization condition reads 



27r' 



ax m — 

1 — X jJL^ 



1. 



(4.14) 



(4.15) 



Within the present approximation (valid only up to 



dx- 



X 



1 

1 — X fi^ 



(4.16) 



In the second term we recognize the familiar expression of wave function correction of the 
state n in old fashioned perturbation theory, namely, Y^'m Ie^^-e\ ■ 
Thus to order a,, 



X 



l+x^ 



1 — X 



(5(1 — x) dy 



Note that ( ^.17| ) can also be written as 

{x,Q' 



X 



8{l 



X] 



—Cfln— 



27r 



1-xU 2 ^ 



X] 



(4.17) 



(4.18) 



which is a more familiar expression. By construction, | <l>2( 

probability den- 
sity. However, this function is singular as x — 1 (gluon longitudinal momentum frac- 
tion approaching zero). To get a finite probability density we have to introduce a cutoff e 
{xgiuon > e), for example. In a physical cross section, this e cannot appear and here we have 
an explicit example of this cancellation. Note that the function P^g = Cjj^^^ + |5(1 — x) 
does not have the probabilistic interpretation since it includes contribution from virtual 
gluon emission. This is immediately transparent from the relation 
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/ 



dxPqq{x) 



0. 



(4.19) 



We also note that the divergence arising from small transverse momentum (the familiar mass 
singularity) cannot be handled properly in the present calculation. This is to be contrasted 
with the calculation of the physical hadron structure function where the mass singularities 
can be properly absorbed into the non-perturbative part of the structure function. 

Let us now explicitly check the longitudinal momentum sum rule for the dressed quark. 
According to the sum rule, 



dx 







2(g) 



» + f: 



2(G) 



X 



,{p I r+(o) I P), = 1. 



(4.20) 



2(P+)2 «^ 

Explicit calculations show that the gluon structure function for the dressed quark target 
^2(0 is given by 

\2 



a i + 
—In—CfX 

271 fX^ ' X 



x] 



From Eqs. (|]T|) and (|]2T|) it follows that 



dxFn(x) 



dx 



[x] + P 



2(G) 



X] 



(4.21) 



(4.22) 



smce 



dx X 



Pgg{x) + Pcgix) 



0. 



(4.23) 



B. Dressed quark structure function from the transverse component 

From BJL expansion and light-front current algebra, it also follows that 
F2{x,Q'') 1 



x 

1 



— I d7]e-''"'Vi 
An J 

rfr7e-%(A;s|V^(r)7V(0)-?^(0)7V(r)|A:s)p. (4.24) 



47rP 



From Eqs. (|4.1| ) and (|4.24|) , it follows that the structure function F2 can be expressed 



not only as a matrix element of the plus component of the bilocal vector current, but also 
the matrix element of the transverse component of the bilocal current. Next we extract 
the structure function F2{x) from the transverse component of the bilocal vector current 
(Eq. ( |4.24]) ). The operator that appears in this equation is 

i^ivh^m = (V^+)t(y)aV-(0) + (^-)t(y)aV(0). (4.25) 

The constrained fermion field ip~ = -r^^a^.^id-^ + gA^) + 'y^m)ip~^. Hence the operator in 
the above equation appears to be higher twist (twist three). Without loss of generahty we 
take the _L direction along the x axis. The structure function can be explicitly written as 
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X 



l_p+ 

8^P 



0^ + 0k^ + 0g\m I P) + h.c, 



(4.26) 



with 



im 



id+ 



a 



id+ 
1 



a 



3^21 



g — [A^ + ^a'A^] 



(4.27) 



First consider contribution from the operator Om- Only potential non-vanishing contri- 
butions are from the diagonal matrix elements for the single quark state and the quark-gluon 
state. Single quark matrix element vanishes because o"^ flips helicity. Diagonal contribution 
from the quark-gluon state also vanishes because of the cancellation between the two terms 
in Eq. ( 4.25 ). Thus the contribution from Om to F2 vanishes. 

Next consider contribution from the operator Of^±. Explicit evaluation leads to 



F2{x,Q' 



X 



+— ^ J dX2 J d^K^ J d'^K2S{l — X — X2)S'^{k^ + K2) 



X \^^2"'\x,Ki;x2,4) 



2 H\+XP^ 



X 



} 



— X] 



+ ^ y dx2 J d^K^ j (i^fi;^(5(l — X — X2) 



since / d k^ki \ ip2 



X 5\Ki + K^) I <|.^^'^^(x,/s:^;X2,4) P } (4.28) 
as a consequence of rotational invariance. Eq. ( |4.28| ) gives the same 



result as Eq. (|4lO|) . 

Lastly we evaluate the contribution from the quark-gluon correlation operator Og. 



F2{x,Q^) _1 
la 2 



X 



2(27r 



,3P1 



E 

o'i,A2 



dy 



(71, A2 



{y,f^^; 1 -y,-i^^) + h.c. 



0. 



(4.29) 



This is because the quark-gluon amplitude $2 has two types of terms: a) terms proportional 
to the quark mass m accompanied by o"-*" which vanish because xt'^'^'Xcr = 0, and b) terms 
proportional to which vanish because of rotational symmetry. Thus the contribution 



from Og to the structure function vanishes. 



From Eq. ( |4.28| ) and Eq. ([4.10| ), it follows that the structure function extracted from Eq. 
( [4.24| ) has the same result given by Eq. ([4.17|) and hence the same parton interpretation as 
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that extracted from Eq. ( [4.1| ). Thus we have exphcitly demonstrated the parton interpreta- 
tion of the transverse component of the bilocal vector matrix element in unpolarized deep 
inelastic scattering. The classification of twist in DIS or other hadronic collision processes 
based on the different components of light-front bilocal operators seems unreliable. 



C. Dressed gluon structure function 

The dressed gluon state can be expanded as 
\PX)g = VAfg{a^{P,X)\0) 



V'2(gq)(-P, A I kiai,k2(J2)h\ki(Ji)d\k2,(y2) \ 0) 
dk-^ d k-^ dJ\i<2 d 



Ai A2 



2(27r)3A;^- ^2{2ny^kt 
ip2{gg){PA I A;iAi,/c2A2)a"^(A;iAi)a^(fc2, A2) | 0)}. 

The factor | is the symmetry factor for identical bosons. 
As before we introduce the boost invariant amplitudes 

VP+ij2{qq){kt , k^) = <i>2iqq){Xi,H^), 
\/P+1p2igg)ikt , k^) = $2(93) 

The qq wave function of the dressed gluon is given by 



^ — — 



x{l—x) 



X 



2(27r) 



a .K 



X 



a .K 



m 



-a — CT" 



— t- 



X 



x(l — x) 



-a 



X- 



The gg wave function of the dressed gluon state given by 

g ^.„„u.x[l-x] i 1 



2(S9) 



2(27r) 



k'Su H 5a + 



-5, 



Hi -I <~>il I :, 

X 1 — X 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



The contribution from the first term in Eq. ( |4.3CI| ) to the gluon structure function for 
the dressed gluon target is given by 



^2(G) 



5{l-x). 



(4.34) 



The contribution to the gluon structure function from the qq component of the dressed gluon 
state is a disconnected contribution which we omit. The contribution to the gluon structure 
function from the gg component of the dressed gluon state is given by 
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r9(3) 



ln^2iV 



1 — X 



1 — X 



X 



X. 



(4.35) 



We define the probability density to a find a gluon with momentum fraction x in the dressed 
gluon, Pgg{x) by 



H h x(l — X 



-1 — X X 

Collecting the three contributions together, we have 

X 



^l(G)(^' = [-^(l - x) + ^ In 9l2N[-^ + ^—^ + x(l - x)]x 



X 



(4.36) 



(4.37) 



The coefficient is determined from the longitudinal momentum sum rule for the 
dressed gluon target, namely, we require. 



/ dxFKx) = dx 
Jo Jo 



1 



2(P+) 



2 9 



Thus we need to evaluate 



2(P+) 



2 9 



{P I ^++(0) I P)g = 1. (4.38) 



(4.39) 



Explicit evaluation leads to 
1 



- (P I 6'++fO) I P) = —In—- I dx 



x^ + (1 -x)2 



(4.40) 



We define the probability density to find a quark with momentum fraction x in a dressed 
gluon as the splitting function Pqcix) 



1 r 



X + (1 — x) 



l + ^\n% [ dxUx' + (1 - xf] + 2N[-^ + 
2tt J ^ 1 — X 



From Eq. ( |4.38| ) we arrive at 

as 

Thus to order a^, we have 



1 — X 



+ x(l — X 



)]4 



(4.41) 



(4.42) 



AT, = 1 - ^ In ^ / rfxf [x^ + (1 - x)^] + 2N[-^ + ^ + x(l - X 
2n J ^ 1 — X X 



)]x}. (4.43) 



Correspondingly, the complete dressed gluon structure function is given by 



Fl^a)i^,Q')=5{l-x) + ^\n 



2-K 

■^^-^ + + x(l - x)]x + ^5(1 - x) 

-XL. X 12 



-5(1 -x)}. (4.44) 
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Including the end point (x — 1) contributions, we define, 
Pgg{x)=2N{' 



1-x 



+ - — - + - x)] + ^(5(1 - x)} - ^6{1 



X 



12 



— x] 



(4.45) 



To the best of our knowledge, this is the first time gluon splitting funtion has been calculated 
using multi-parton wave-funtions. There exist some discussions in the literature regarding 
the calculation of DIS splitting funtions using the language of multi-parton wave-funtions 
mainly due to Lepage and Brodsky [||. But for the gluon splitting funtion, they have simply 
quoted the result from Altarelli-Parisi paper It is easy to verify that 



dx X 



0. 



(4.46) 



V. POLARIZED DRESSED PARTON STRUCTURE FUNCTIONS 

In this section we discuss in detail the calculations of polarized structure functions of a 
dressed quark target in light-front perturbation theory. 



A. Chirality structure function 

The chirality structure function of a dressed quark is given by 

= ^ / rfr7e-%(A;s|V^(r)7+75V^(0)+V^(0)7+75^(r)|fcs)p. 



(5.1) 
(5.2) 



On the light-front, chirality and intrinsic helicity of a fermion coincide. A direct calculation 
of gi{x) from Eq. ( |5.2|) in the free quark helicity state immediately leads to the well-known 
solution: 



9iix) = -^5(1 -x). 
By a calculation similar to that of F2, for gi, we have for the dressed quark. 



9i{x,Q') = '-^l6{l~x) + ^Cjln% 



27r 



(5.3) 



(5.4) 



which has the same form as Eq. ( 4.18 ). Therefore, the splitting function for gi is the same 
as that for F2. 
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B. Transverse polarized structure function 



The transversely polarized structure function of a dressed quark is given by 

9t{x,Q') = ^ fd^e-'^'^A, (5.5) 

= ^J dve-''^%{ksmn{l' - ^l^)l^m + h.c.\ks),. (5.6) 

For convenience, we take the polarization along the x-direction. The transverse polarized 
quark target in the x-direction can be expressed in terms of helicity states by 

k^, S') = k^, ]) ± k^, j)) (5.7) 

with = ±m^, and is the renormalized quark mass. The operator (see Eq. (|5.6| )) 

?^(r)7±75V'(0) = + O''^ + (5.8) 

where 

O^ = #Ur)(4±(r)^7±-7±^/±(0))75^+(0) (5.9) 

where m and g are the quark mass and quark-gluon coupling constant in QCD, and A_[_ = 
A"^Ta is the transverse gauge field. 

Without the QCD correction (i.e., for the free quark state), it is easy to show that 



grix) = gl^ix) = |^5(1 - x) = |<5(1 - x), g^^-{x) = = ^|^(x). (5.10) 



Here rriq/S^ = 1 since the renormalized mass is the same as the bare mass at the tree level 
of QCD. We see that only the quark mass term contributes to gT in Eq. ( |5.9|) . The quark 
transverse momentum term alone does not contribute to gx since it cannot cause helicity 
flip in the free theory. This result indicates that physically the dominant contributions to 
gT is not controlled by the twist classification. Thus, for the free quark, we have 

92ix) = grix) - gi{x) = 0. (5.11) 

It is obvious that for free theory, the Burkhardt-Cottingham (BC) sum rule is trivially 
obeyed. But as we can see (as it has also been previously noticed) the Wandzura-Wilczek 
relation |18|, 



92ix) = -g,ix) + I'dy^, (5.12) 

y 



16 



is not satisfied in free tlieory. 

Next, we consider tlie QCD corrections up to order a^, where the quark-gluon interaction 
is explicitly included. We find that all the three terms in Eq. ( |5.9| ) have nonzero contribution 
to Qt, 



el rUg 



2tt 



1 — X 



— 6(1 — x) [ dx — i I 

^ 'Jo 1 - X' J J 



g^^{x,Q'^) 
g^{x,Q'^) 



el nig as , ^ 
Q'6il-x) 



(^q^gas 



(5.13) 
(5.14) 
(5.15) 



from the dressed quark wave function normalization constant A/" in Eq. ( [4.4|) (corresponds 
to the virtual contribution in the standard Feynman diagrammatic approach). It shows that 
up to the order as, the matrix elements from Ok^^ (quark transverse momentum effect) and 
Og (quark-gluon interaction effect) in Eq. ( |5.9| ) are also proportional to quark mass. In other 
words, the transverse quark momentum and quark-gluon coupling contributions to grix, Q^) 
arise from quark mass effect. Explicitly, these contributions arise from the interference of 
the rriq term with the non-mg dependent terms in the wave function of Eq. ([4.9| ) through 
the quark transverse momentum operator and the quark-gluon coupling operator in the gx 
expression. This result is not surprising since, as we have pointed out, the pure transverse 
polarized structure function measures the dynamical effect of chiral symmetry breaking [jl9| . 



Physically only those interferences related to quark mass can result in the helicity flip (i.e., 
chiral symmetry breaking) in pQCD so that they can contribute to gT{x,Q^). From this 
result, we may see that only the operators themselves or their twist structures may not give 
us useful information about their importance in the determination of structure functions. 
Combining the results together, we obtain 



9t{x,Q^ 



x] 



a 



27r 



1 + 2a; 



x] 



25(1 



X 



(5.16) 



Note that in the above solution, m„ is the bare quark mass, while the dressed quark polar- 



ization 5"^ 



rri' 



and up to order ag 



1 + —a^Cf m 



47r 



(5.17) 



We must emphasize that on the light-front there are two mass scales in the QCD Hamil- 
tonian, one is proportional to which does not violate chiral symmetry, and the other is 
proportional to nig which we discuss here and is associated with explicit chiral symmetry 
breaking in QCD [13]. An important feature of light-front QCD is that the above two mass 
scales are renormalized in different ways even in the perturbative region. The renormaliza- 
tion of ml in pQCD is different from the above result, the details of which can be found in 
our previous work B. With this consideration, we have 
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gT{x,Q^ 



l + 2x 



X 



'l-x) 



+ 2^(1-^) 



(5.18) 



The final result is independent quark mass but we have to emphasize again that one must 
start with massive quark theory [Q. Otherwise, there is no definition for qt at the beginning. 
This result has a close analogy in helicity fiip process at high energy in Quantum Electro 
Dynamics ||2l|]. Cross section for such process, which vanish in the chiral limit according to 
naive arguments, indeed is non-vanishing if one lets the electron mass to go zero at the end 
of the calculation as Lee and Nauenberg pointed out long time ago p2 |. 
Thus, up to of the oder as, we find g2 for a quark target 



92ix, Q') = -^-^Cf \n^2x- 5(1 - x 
It is easy to check that the above result of g2{x, Q"^) obeys the BC sum rule. 



/ dxg2{x,Q'^) = 0, 
Jo 



(5.19) 



(5.20) 



as is expected. Our results also show that Wandzura-Wilczek relation is strongly violated 
in perturbative QCD. 



VI. STRUCTURE FUNCTION OF HADRON: PARTON PICTURE, SCALE 
EVOLUTION AND FACTORIZATION 

As we have schematically discussed in Sec. II (also see 0), the nonperturbative con- 
tribution to the structure functions and the scaling violations from the perturbative QCD 
corrections can be unified and treated in the same framework in our formalism. In this 
section, we shall address the issues associated with scaling violations in the structure func- 
tion of the "meson-like" bound state and explicitly demonstrate the validity of factorization 
outlined in Sec. II. 



A. Parton picture 

Let us first discuss the emergence of parton picutre for the structure function of a com- 
posite state. We expand the state | P) for qq bound state in terms of the Fock components 
qq, qqg, ... as follows. 

di k-^ r dk'2 d k<2 



ai,<72^ .J2(27r)3A;^ ^ J2{2nYk^ 



MP I ki,cTi; A;2,a2)V2((27r)3p+(5='(P - k, - k2)b\k,,a,)d\k2,a2) \ 0) 
^ -^ I' dk-^ d k-^ ^ dk<2 d k<2 ^ dk^ d k^ 



a^,a,[xj ^2{2'Kfkt ^ ^2{2-Kfkt ^ ^2(2^0^ 



i;^{P I A;i, di; k^, a^, k^, \^)p{2'ixYP^8\P -k^-k^- k^) 

b\ki,ai)d\k2,a2)a\ks,Xs) \ 0) 
+ ... . (6.1) 
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Here V'2 is the probability amplitude to find a quark and an antiquark in the meson, i/j^ is 
the probability amplitude to find a quark, antiquark and a gluon in the meson etc. 

As in Sec. IV we evaluate the expression in Eq. ( ^4.1| ) explicitly. The contribution from 
the first term (from the quark), in terms of the amplitudes 

P^^,{kt,ki) = ^,{x,,4), (6.2) 



and so on, is 
F|(x) 



X 



/ dX2 / d^K-^ I (i^/t^(5(l — X — X2)(5^(Ki + K2) I $2^''^^(x, K^; X2K^) P 
0-1,0-2 

+ ^ y dx2 j dx^ j d^K^ j d^K2 j d^K^5{l — x — X2 — x^)5'^{ki + K2 + K3) 



01,0-2, A3 ' 



$^-'^^''^(x,«^;x2,4;x3,4) r+- • (6.3) 



Again, the partonic interpretation of the F2 structure function is manifest in this expression. 
Using different techniques and approximations, the same result has been also obtained by 
Brodsky and Lepage 0. 

Contributions to the structure function from the second term in Eq. ( [4.1|) is 



[X 



X 



/ dX2 / (i^K^ I (i^K^5(l — X — X2)5^(ki + K2) I $2^''^^(X2, K^; X, K^) 1^ 

01,0-2"' 

+ ^ y dX2 j dXs J d'^K^ J d'^K2 J d'^K^5{l — X — X2 — X-i)5'^{Ki + K2 + K3) 
01 ,0-2, A3 

I n"''"''{x2,4;x,^^t,x3,4) r +••• • (6.4) 

The normalization condition guarantees that 

. f?W+^lM|=2 (6.5) 

J X X 

which refiects the fact that there are two valence particles in the meson. Since the bilo- 
cal current component ^ involves only fermions explicitly, we appear to have missed the 
contributions from the gluon constituents altogether. Gluonic contribution to the struc- 
ture function F2 is most easily calculated by studying the hadron expectation value of the 
conserved longitudinal momentum operator P"*". 

From the normalization condition, it is clear that the valence distribution receives contri- 
bution from the amplitudes $2, ^3, ■■■ at any scale /i. This has interesting phenomenological 
implications. In the model for the meson with only a quark- antiquark pair of equal mass, 
the valence distribution function will peak at x = |. If there are more than just the two 
particles in the system, the resulting valence distribution will no longer be symmetric about 
X = 2 as a simple consequence of longitudinal momentum conservation. 

The equation ( |6.3D as it stands is useful only when the bound state solution in QCD 
is known in terms of the multi-parton wave functions. The wave functions, as they stand, 
span both the perturbative and non-perturbative sectors of the theory. Great progress in 
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the understanding of QCD in the high energy sector is made in the past by separating 
the soft (non-perturbative) and hard (perturbative) regions of QCD via the machinery of 
factorization. It is of interest to see under what circumstances a factorization occurs in the 
formal result of Eq. (|6.3|) and a perturbative picture of scaling violations emerges finally. 
We shall explicitly address this issue in the following section. 



B. Perturbative picture of scaling violations in a bound state 

To address the issue of scaling violations in the structure function of the "meson-like" 
bound state, it is convenient to separate the momentum space into low-energy and high- 
energy sectors. Such a separation has been introduced in the past in the study of renormal- 
ization of bound state equations p3 in light-front field theory. The two sectors are formally 



defined by introducing cutoff factors in the momentum space integrals. How to cutoff the 
momentum integrals in a sensible and convenient way in light-front theory is a subject under 
active research at the present time. Complications arise because of the possibility of large 
energy divergences from both small and large k-^ regions. In the following we investigate 
only the effects of logarithmic divergences arising from large transverse momenta, ignore 
subtleties arising from both small x {x —>■ 0) and large x {x —>■ 1) regions and subsequently 
use simple transverse momentum cutoff. For complications arising from x — 1 region see 
Ref. i. 



1. Scale separation 

We define the soft region to be < /i and the hard region to be < k"*- < A, 
where n serves as a factorization scale which separates soft and hard regions. Since it is an 
intermediate scale introduced artificially purely for convenience, physical structure function 
should be independent of /i. The multi-parton amplitude $2 is a function of a single relative 
transverse momentum and we define 

r $i, < < /i. 

The amplitude $3 is a function of two relative momenta, Ki and K2 and we define 

r < < /i 

^ I $f , < ni < fi, < 4 < A 

^ I fi< < A, < < fx ^ ^ ^ 

[<^>t^ fi<Ki,4<A 

Let us consider the quark distribution function q(x) = -^^^ defined in Eq. (|6.3|) . In presence 
of the ultraviolet cutoff A, q{x) depends on A and schematically we have, 

g(^,A^)=Er*2+Err*3- (6.8) 

Jo Jo Jo 

For convenience, we write. 
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q{x, A^) = q^ix, A^) + g3(x, A^). (6.9) 

where the subscripts 2 and 3 denotes the two-particle and three-particle contributions re- 
spectively. Thus, schematically we have, 



rA 



A 

h |2 
2 



r^l rA 






2 


Jo Jfi 




/•A rA 








J fj, J ^ 





r-A rfl 

I I rf)hs |2 



Jfi Jo 

+ E/ / (6-10) 

We investigate the contributions from the amplitudes $3'* and ^'^^ to order in the follow- 
ing. 



2. Dressing with one gluon 



We substitute the Fock expansion Eq. ( |6.1|) in Eq. ([4.3|) and make projection with a 
three particle state b\ki, ai)S{k2, £12)0^(^3, as) \ 0) from the left. In terms of the amplitudes 
$2, $3, we get, 

$^^"2^^(a;, Ki; X2, k^.I-x- X2, K3) = Mi + M2, (6.11) 
where the amplitudes 

Ml = h-) JL — T° ^ $^^^^(1 - X2, -4; ^2, 4), (6.12) 

^ V2(27r)3 VI-X-X2 

and 

M2 = ^^==T" ^ F2 4; 1 - X, -4) (6.13) 

^ y'2(27r)3 VI - a; - X2 



with the energy denominator 

,2 _ 

x X2 1 — X — X2 



E^lM^--!!l±Ml-!!l±Ml (6.14) 



and the vertices 

^1 — X — X2 a; 1 — X2 

and 

= E \-r^ — ^Ast^i^ ^ (6.16) 

,1 — a; — X2 a;2 1 — x 21 

■^2 
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3. Perturhative analysis 



For hard and ^2 soft, + Kg ~ Hi and the multiple transverse momentum integral 
over $3 factorises into two independent integrals and the longitudinal momentum fraction 
integrals become convolutions. The contribution from Ali to $3 is, 



^3,1 ^{x,K^;x2,K2;l-x- X2, -K2 ) = , r 



^2 



2(27r)3 1 - a:^2 (k^)^ 



a; — a;2 a; 

$2""'(i-^2,-4;^2,4)- (6-17) 

Thus the contribution from A^i to the structure function is 

A2 r^dy 



(6.18) 



where 




l + (-)2 

1 - ^ 



(6.19) 



For the configuration k^ hard, k^ soft, contribution from AI2 does not factorize and the 
asymptotic behavior of the integrand critically depends on the asymptotic behavior of the 
two-particle wave function $2- To determine this behavior, we have to analyze the bound 
state equation which shows that for large transverse momentum $2('t"'") ~ (^^rya- Thus 
contribution from M.2 for scale evolution is suppressed by the bound state wave function. 
Analysis of the interference terms (between }A\ and AI2) shows that their contribution also 
is suppressed by the bound state wave function. 

For the configuration k\; soft, k^ hard, contributions from M.\ and the interference terms 
are suppressed by the wave function. Contribution from M.2 factorises both in transverse 
and longitudinal space and generate a pure wave function renormalization contribution: 

(6.20) 

We have seen that even though the multi-parton contributions to the structure function 
involve both coherent and incoherent phenomena, in the hard region coherent effects are 
suppressed by the wave function. 



^. Corrections from normalization condition 

In the dressed quark calculation, we have seen that the singularity that arises as x' ^ 1 
from real gluon emission is canceled by the correction from the normalization of the state 
(virtual gluon emission contribution from wave function renormalization). In the meson 
bound state calculation, so far we have studied the effects of a hard real gluon emission. In 
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this section wc study the corrections arising from the normahzation condition of the quark 
distribution in the composite bound state. 

Collecting all the terms arising from the hard gluon emission contributing to the quark 
distribution function, we have, 

+ ^Cfln^q,{x, Pi') J dyP{y). (6.21) 

We have a similar expression for the antiquark distribution function. 

The normalization condition on the quark distribution function should be such that there 
is one valence quark in the bound state at any scale Q. We choose the factorization scale 
fjL — Qo. Let us first set the scale A — Qo. Then we have (in the truncated Fock space) 

/' dx q2{x, Ql) + r dx qsix, Ql) = 1. (6.22) 
Jo Jo 

Next set the scale A = Q. We still require 

r dx q2ix, Q') + C dx qsix, Q^) = 1. (6.23) 
Jo Jo 

We note that the evolution of gs requires an extra hard gluon which is not available in the 
truncated Fock space. Thus in the present approximation q3{x,Q') — qz{x,Q^)- 
Carrying out the integration explicitly, we arrive at 

dx q2{x, Ql) [1 + ^Cfln^^ J dyP{y)] + £ dx q^{x, Q^) = 1 (6.24) 

Thus we face the necessity to "renormalize" our quark distribution function. Let us 
define a renormalized quark distribution function 

Q^o) = ?2(x, Q^)[l + 2-^C//n^ / dyP{y)\ (6.25) 

so that, to order a^, 

C dx q^{x, g^) + dx qsix, Ql) = L (6.26) 
Jo Jo 

We have, 

g2(x, Ql) = q^{x, Ql){l - 2^^^/^^^ I dyP{y)]. (6.27) 
Collecting all the terms, to order as, we have the normalized quark distribution function. 
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q{x,Q'')=q^{x,Ql) 

+ ^Cfln% f dy q^{y, Ql) f dz 5{zy - x) P{z) 
Zn C^Q JO JO 

+ g3(x,g') (6.28) 

with P{z) = P{z) - 6{z - 1) dyP{y). 

We see that just as in the dressed quark case, the singularity arising as x ^ 1 from 
real gluon emission is canceled in the quark distribution function once the normalization 
condition is properly taken in to account. From this derivation we begin to recognize the 
emergence of the Altarelli-Parisi evolution equation. 

C. summary 

In this section we have carried out an analysis of the scale evolution of structure functions 
of a meson-like composite system. We have separated the parton transverse momenta into 
soft and hard parts. The three body wave function which is a function of two relative 
momenta has soft, hard and mixed components. The mixed components of the three body 
wave function which are functions of soft and hard momenta are responsible for the scale 
evolution of the soft part of the structure function in perturbation theory. 

In the analysis with wave functions, there are two contributions to the three body wave 
function: One where the gluon is absorbed by the quark and second where the gluon is 
absorbed by the anti-quark (spectator). There appears a non- vanishing contribution when 
the hard gluon is absorbed by the anti-quark. This corresponds to the transition caused 
by the interaction Hamiltonian when the active parton remains soft, while a hard spectator 
makes transition to a soft spectator state. This leads to wave function renormalization of 
the spectator anti-quark but this is eventually canceled by the normalization condition as 
discussed in detail in Sec. VIB.4. This justifies a posteriori the prescription given in Sec. II 
that we need to keep only those terms in P"(^) which cause transitions involving the active 
parton. 

In the wave function analysis, there are also contributions that are omitted a priori in 
the calculational scheme which lead to factorization in Sec. II. All of these contributions 
are suppressed by the asymptotic behavior of the bound state wave function as we have 
explicitly shown. In summary, the detailed analysis carried out with the help of multi- 
parton wave functions in Sec. VIB justifies the approximations made in Sec. II which lead 
to the emergence of factorization to all orders in perturbation theory and to the simple scale 
evolution picture. 

VII. CONCLUSION 

We have shown that a perturbative analysis in the light-front Hamiltonian formalism 
leads to the factorization scheme proposed in Ref . . It is shown that the scaling violations 
due to perturbative QCD corrections can be rather easily addressed in this framework by 
simply replacing the hadron target by dressed parton target and then carrying out a sys- 
tematic expansion in the coupling constant a.^ based on the perturbative QCD expansion of 
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the dressed parton target. The calculational procedure utihzes techniques of old-fashioned 
perturbation theory main ingredients of which are transition matrix elements and energy 
denominators. 

The main advantage of the present method can be summarized as follows. The bilocal 
currents are defined in the light-front gauge A'^ = and since the bilocality is only in the 
light-front longitudinal {x~) direction, the path-ordered exponential between fermion field 
operators in the bilocal current is replaced by unity [irrespective of which component of 
the current is considered). This results in an extremely simplified operator structure and 
a straightforward parton picture. Further, the calculations do not employ Feynman prop- 
agators and as a result we encounter neither the usual problems associated with using a 
non-covariant gauge in a covariant calculation nor the problems associated with the un- 
physical pole of the propagator. The calculations are straightforward and 75 or presence of 
quark masses pose no special problem. The physical picture is very clear at every stage of 
the calculation. Also the regularization scheme used in this framework for perturbative con- 
tribution can be directly applied to the construction of hadronic bound states which is the 



major topic of current research on light-front field theory p^JT^ . Thus, once the light-front 
bound state structure are found, a complete theoretical understanding of structure functions 
can become possible. 

In addition, the approach uses probability amplitudes rather than probability densities 
and hence interference effects are easy to handle. Exploiting this feature, we have clarified 
the parton interpretation of the matrix elements of the transverse component of biocal vector 
and axial vector currents. We have presented real and virtual corrections to the structure 
functions F2 and gi for a dressed quark and gluon in a transparent manner. The splitting 
functions are extracted and the longitudinal momentum and helicity sum rules are verified 
explicitly to order a^. 

We have carried out, with the help of multi-parton wave functions, a detailed analysis 
of the scale evolution of the structure function of a composite system which justifies the 
approximations made in Sec. II which lead to the emergence of factorization to all orders 
in perturbation theory and to the simple scale evolution picture. A complete fourth order 
calculation is necessary to establish the viability of the new approach for the perturbative 
domain. Such a calculation is presently under way. Also, the main contribution to DIS 
structure functions, the nonperturbative QCD dynamics, is also in progress. We shall leave 
the discussion of these topics for future publications. 
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